Goal:

In this experiment I am observing a co-evolving population (survival is determined by a random draw) with a non-interactible trait. I expect to see no spatial phenotype correlation

Questions:

What happens when a co-evolving trait is not interitable?

Background

## All cor, lit, and grid files exist!
## This program will now end!

Important file names

file_name= “Snake_1e-11_5_Newt_1e-11_5_0” file_name= “Snake_1e-08_0.005_Newt_1e-08_0.005_3” inter_tall_tall_file inter_cor_tall_file inter_inter_grid_tall_file

Mean Phenotype Whole Simulation

notes

Mean Phenotype Whole Simulation

notes

Phenotype differences

Table of average Differences

##                    Group.1           x
## 1  1e-08_0.005_1e-08_0.005 -0.18729482
## 2   1e-08_0.005_1e-09_0.05  0.13976018
## 3    1e-08_0.005_1e-10_0.5  0.50759825
## 4      1e-08_0.005_1e-11_5  1.07418888
## 5   1e-09_0.05_1e-08_0.005 -0.11800217
## 6    1e-09_0.05_1e-09_0.05 -0.08351749
## 7     1e-09_0.05_1e-10_0.5  0.57878556
## 8       1e-09_0.05_1e-11_5  0.74688826
## 9    1e-10_0.5_1e-08_0.005 -0.61776623
## 10    1e-10_0.5_1e-09_0.05 -0.59585200
## 11     1e-10_0.5_1e-10_0.5  0.02503837
## 12       1e-10_0.5_1e-11_5  0.45853498
## 13     1e-11_5_1e-08_0.005 -0.97525745
## 14      1e-11_5_1e-09_0.05 -0.47443100
## 15       1e-11_5_1e-10_0.5 -0.42495875
## 16         1e-11_5_1e-11_5  0.04003739

Explain

Connection between higher phenotype and population

Notes

Phenotype differences

Phenotype & Populationsize differences

explain

Correlation

notes

After looking at this figure, I can see that there is a larger range of correlation values. The can sometimes be positive, near 0, or negative. Very few reach (or come near) the real newt-snake correlation. How do you find the correct time slice to look for spatial correlation? How would these phenotypes become spatially correlated and then become not spatially correlated.

Correlation Histograms

In order to understand how spatial correlations where changing with time I took 5,000 generation time slices to look at all four trials correlation values. Each color is a different trial per GA combination. The histogram values are stacked.

Plot 1

Plot 2

Plot 3

Plot 4

Plot 5

Plot 6

Plot 7

Plot 8

Plot 9

Plot 10

exp

Correlation across time

notes

Random 1

Random 2

## [1] "pattern 1e-08_0.005_1e-09_0.05_1"
## [1] "Cor between average snake pheno and local cor -0.177345827911184"
## [1] "Cor between average newt pheno and local cor -0.177394682818902"
## [1] "Cor between average dif pheno and local cor 0.118831497227703"
## [1] "Cor between newt pheno and snake -0.270520115883644"

Random 3

## [1] "pattern 1e-09_0.05_1e-11_5_0"
## [1] "Cor between average snake pheno and local cor 0.16777243970841"
## [1] "Cor between average newt pheno and local cor -0.00726983202923471"
## [1] "Cor between average dif pheno and local cor 0.133559933250909"
## [1] "Cor between newt pheno and snake 0.600181538612102"

I still wonder why the correlation values go up and down, because the mean newt and snake phenotype keeps going up. I wonder if the individuals might be moving around too much. I wonder what would happen on a larger map?

What happens over time (looking at the beginning, middle, and late part of my simulations)

This next section is just getting a glimpse at how newt & snake phenotype and population size differ over time. The populations start off with about 250 individuals each. Each individual has a different genetic background created from msprime.

Pheno Beginning

Pheno Middle

Pheno End

Dif Beginning

Dif Middle

Dif End

explane

Summary

notes

Early-Sim Population Size Summary

Mid-Sim Population Size Summary

Late-Sim Population Size Summary

Early Difference Summary

Mid Difference Summary

Late Difference Summary

By Snake GA (Early)

By Snake GA (Mid)

By Snake GA (Late)

By Newt GA (Early)

By Newt GA (Mid)

By Newt GA (Late)

expla

Heatmap

notes

Population Size (Early)

Population Size (Mid)

Population Size (Late)

Phenotype (Early)

Phenotype (Mid)

Phenotype (Late)

exp

What is up with the correlations

notes

Early Simulation Correlation

Mean

## [1] 0.08941538

SD

## [1] -0.05684025

Max Phenotype

## [1] -0.05620503

Min Phenotype

## [1] -0.1811683

Population Size

## [1] -0.2839627

Middle Simulation Correlation

Mean

## [1] -0.2447956

SD

## [1] 0.2091648

Max Phenotype

## [1] -0.19622

Min Phenotype

## [1] 0.01225089

Population Size

## [1] 0.6059479

Late Simulation Correlation

Mean

## [1] 0.07210036

SD

## [1] 0.4470126

Max Phenotype

## [1] 0.241293

Min Phenotype

## [1] 0.08881336

Population Size

## [1] 0.440485